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Abstract 

The SLAC data on the p, d and n polarized structure functions are fairly well reproduced with 
and without the contribution of the anomaly. The results are compared with a previous study 
based mainly on SMC data. The implications on the solution of the spin-crisis are discussed. 



Wn = ^ + ^)(l-^)=^(l-^), (3) 



The EMC experiment Q on polarized deep inelastic scattering muon-proton at <Q 2 >= 10.7 
GeV 2 led to the result 

I = ( l g p 1 {x)dx = 0.126 ±0.010 ±0.015 (EMC/SLAC), (1) 
Jo 

smaller than the value predicted by the Ellis and Jaffe sum rule Q 

F D 

= 0.185 ±0.001, 2) 

2 18 K ' 

and gave rise to the "spin crisis" . An explanation of this result has been proposed in terms of the 

anomaly of the axial vector current related to the gluon polarization AG Q, but, to recover the 

EMC result, the rather high value of AG ~ 5 at <Q 2 >= WGeV 2 was necessary, as pointed out 

in Ref. @, to be compensated by a large negative orbital angular momentum <L Z >. This large 

component of the proton spin carried by the gluon and this compensating mechanism is far from 

being natural, but the Bjorken sum rule || including first order QCD corrections, 

would be obeyed and the generally accepted framework to describe deep inelastic phenomena would 
not be spoiled. 

The EMC result, assuming the validity of the Bjorken sum rule, would imply a large negative 
value for the neutron Ellis and Jaffe sum rule 

In= [ gUx)dx, (4) 
J 

while the measured value at SLAC by E142 at <Q 2 >= 2GeV 2 || with a polarized He 3 target, is 

I n = -0.022 ±0.011, (5) 

in good agreement with the prediction of Ellis and Jaffe in the absence of the anomaly contribution: 

— — -D = —0.022 ± 0.004. (6) 
3 9 y J 

Higher twist operators have been advocated |7| to play an important role in the interpretation 
of the E142 data, which has been taken at the rather small <Q 2 >= 2GeV 2 . 

More recently, at CERN |8| and SLAC Q new deep inelastic scattering experiments have been 
performed with polarized proton and deuteron targets. An analysis of the available data has been 
achieved by Gehrmann and Stirling 0] to describe the gi(x) distributions. In the framework of 
the explanation of the defect in the Ellis and Jaffe sum rule for I p by means of the gluon anomaly, 
they try to describe g\{x) and g\{x) in terms of the contributions of the valence quarks, Au v (x) 
and Ad v (x), and of AG(x): 

gl{x,Q 2 ) = lAu v (x,Q 2 ) + ^-Ad v (x,Q 2 )-\^^AG(x,Q 2 ), (7) 
y 18 6 In 

9i{x,Q 2 ) = ^Au v (x,Q 2 ) + ^Ad v (x,Q 2 )-^^^-AG(x,Q 2 ), (8) 



2 



with 

xAu v = r] u A u x au (1 - x) K (1 + 7n x) 
xAd v = r] d A d x ad (l-x) bd (l+j d x) 
xAG = r] G A G x aa (I - x) b ° (I + -f G x) 
and A q = A g (a q , b q ,^ q ) (q = u,d,G) in such a way that 



(9) 



/ Aq(x)d. 
Jo 



x = rj q . 



(10) 



At Qq = 4:GeV 2 , lower limit for the data considered in order to avoid the uncertainty on the 
higher twist contributions, the rj's are fixed by [O (a s (4GeV 2 ) = 0.2879): 



Vu(Ql) 



2F{Ql) = 2 {1 

F{Ql) - D(Ql) : 
-0.294 ±0.011 



3a s (Q§) 



5-7T 



p 2a s {Ql) D 



15vr 



0.848 ±0.016 



(11) 



5ir 



F 



na s (Q 2 ) 



15vr 



D 



(12) 



rj G (Ql) = 1.971 ±0.929. (13) 

With theoretical and phenomenological motivations they require 

a u = a d a G = 1 

b d = 6 U ±1 = 4.64 (14) 

1u = Id 

and they obtain a very good fit to the data, x 2 /Ndf = 0.63, corresponding to the parameters 

a u = 0.46 ± 0.15, 

j u = 18.36 ± 14.49, (15) 
b G = 7.44 ± 3.52, 

with the choice 7g = 0. 

As the authors stress, the bulk of the SLAC data falls outside their fit region, excluding almost 
all the E142 neutron data. Given the large uncertainty on the SMC deuteron data, one can conclude 
that the constraint coming from g\{x) is not very strong and the fair description obtained for gi(x) 
is not surprising, once r] u , r] d , and r\ G are chosen consistently with the experimental value found for 
I P . 

In order to test the Bjorken sum rule, which is the most appealing target of experiments on 
polarized structure functions, the value of /„ should be measured and one should try to reproduce 
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Concerning the exclusion of data below Q 2 = 4 GeV 2 , because of the poor theoretical knowledge 
on the higher twist contributions, it is worth observing that the values of I p measured at <Q 2 >= 
3 GeV 2 at SLAC §, 0.127 ± 0.004 ± 0.010, and at <Q 2 >= WGeV 2 at CERN |], 0.136 ± 0.011 ± 
0.011, are consistent, showing that for the proton one may safely neglect higher twist contributions. 

Here we try to describe the SLAC data on proton and deuteron at <Q 2 >= 3 GeV 2 and on 
neutron at <Q 2 >= 2GeV 2 in terms of Au v , Ad v , and AG, given by Eq. (||), with the values of 



rj u and % scaled down to Qq = 3 GeV 2 . We are using the values of F and D of Ref . 1 11 ] and 
a s (3GeV 2 ) = 0.35 ± 0.05, 

fj u = 2F = 0.835 ± 0.022 , (16) 
fj d = F-D = -0.276 ±0.020, (17) 

and 

6vr /2 1 \ 

to recover the measured value for I p = 0.127 ± 0.004 ± 0.010 ||. 

As in Ref. |[(]] we fix a u = ad, J u = Id, clq = 1, and further, we require 

b u > 1, b d > 3, b G > 5, lu = ld > 0, and 7G > 0. (19) 

The lower limits for b u and b^ are smaller than the values chosen for these parameters in Ref. 
[|10[| , as shown in our equation (p^): the smaller one for b u is taken since v} dominates at high x. 
The lower limit on be of the order of the corresponding exponent for G in Eq. (15), has also the 



purpose to prevent AG(x) from carrying too much proton momentum (the momentum carried by 
AG should of course be less than the momentum carried by G). In order to have a definite sign 
for each spin asymmetry we do not allow values of 7 < —1. In the range (—1,0) there is a strong 
correlation between the values of 7 and b (in particular 7 = — 1 and b give the same function than 



7 = and 6+1). The large error in the determination of j u in Ref. [10] shows that this parameter 
does not play a crucial role and therefore we assume all the j's to be positive. 
To keep into account the D-state admixture for deuterium, we take 

gt(x) = \(l-\u JD \(gl(x)+g n 1 (x)) , (20) 



where lo d = 0.058 |2|. 



The best fit to the SLAC polarized structure functions of the theoretical expressions given by 
Eqs. - ([lOD and (|i~6|)-(20) is found, with a x 2 /Ndf = 1-8, corresponding to the following values 
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for the parameters {clg = 1): 



a u = ad 


= 0.57 


± 


0.03 


K 


= 2.0 


± 


0.3 


bd 


= 3.0 


± 


0.1 


bo 


= 20 


± 


1 


lu = Id 


= 1.0 


± 


0.8 


1G 


= 0.0 


± 


1.0 



(21) 



The resulting predictions are compared with SLAC data in Figs. 1, 2, 3, showing a tendency to 
get large negative values at small x for gf, gf, just in the region where the negative contributions 
of AG are more important. The behaviour of xAu(x), xAd(x), and xAG(x) is described in Fig. 4. 
If one leaves t\q free, one gets a better fit, x 2 /Ndf = 1-0, with (clg = 1) 



a u — ad 




0.40 


± 


0.08 


K 




1.8 


± 


0.2 


b d 




3.0 


± 


0.5 


bo 




20 


± 


15 


lu = Id 




5.3 




+5.4 
-2.7 


1G 







± 


18 


t)g/vg 




0.48 


± 


0.09 


d rj d free, 


the best 


fit, with x 2 /Ndf 


a u = a d 




0.98 


± 


0.09 


K 




1.8 




+0.5 
-0.2 


b d 




4.9 


± 


0.8 


lu = Id 







± 


3 


Vu/fju 




0.76 


± 


0.03 


ml 'm 




0.93 


± 


0.08 



(22) 



(23) 



The predictions corresponding to the values in Eqs. ( |22|) and (P3j) are also compared with the 
SLAC data in Figs. 1, 2, 3. The parameters given in Eqs. ( |2~I|) and ( p3|) are reported in Table 



[y together with the parameters of reference [10], hereafter referred as Fit A, Fit B, and Fit C 
respectively. It is worth stressing that b u and bd in Fit A come out smaller than the values fixed in 
Ref. [10] , while be comes out larger than the corresponding value found in Ref. [10] and the total 
proton momentum carried by AG in Fit A is about twice smaller than in Fit C. 

The values of rj u and r/d of Fit B support a different interpretation for the defect in the Ellis and 
Jaffe sum rule for the proton, which has been related to the defect in the Gottfried sum rule and 
to the role of Pauli principle for parton distributions [O]. It amounts to say that the v) parton, 
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largely the most abundant in the proton, receives less contribution from the sea because its energy 
levels are almost completely occupied. Indeed, the role of Pauli principle in parton distributions 
explains, in a natural way, the dominance of iJ at high x, shown by the increase towards 1 of the 
asymmetry parameter A^{x) at high x and also the fast decrease in the same limit of the ratio 
F^ixj/F^ix). This was known since the time of the phenomenological analysis performed by Field 
and Feynman fl4}| , who observed that Pauli principle may be responsible for less uu pairs than 
dd in the proton, just due to the presence of two valence u quarks and only one d quark. The 
defect in the Gottfried sum rule confirms their conjecture and it is natural to assume that the v) 
quark rather than the quark receives less contribution from the sea. The fact that rj u , reported 
in Eq. (|23|), is smaller by 25% than the expected value, while the reduction of % is less relevant, 
complies with this picture. 

The importance of the role of Pauli principle in the parton distributions therefore suggests to 
describe them as Fermi-Dirac functions [15, |lG| ]: 



p(x) = f{x) 



. x — x(p) . 
exp i } + 1 



n -1 



(24) 



where x, x(p), and f(x) play the role of the "temperature", "thermodynamical potential", depend- 
ing on the flavour and on the spin of the parton, and "weight function" for the density of levels in 
the x variable. 

Two models have been proposed with the common feature of relating the shape of the distri- 
butions of the different partons to their first moments (broader shapes for larger first moments). 
In one of them the polarized distributions are found, without additional parameters, from the 
unpolarized distributions with the assumptions: 

Au(x) + Au(x) = u(x) + u(x) — d(x) — d(x) (25) 
Ad(x) = (F-D) (d(x) - d(x)) (26) 
Ad(x) = (27) 

In the second one the potentials of u and d of both helicities are put as parameters in order 
to describe unpolarized and polarized available distributions. 

Eq. (|^) would imply, with only u, d and their antiparticles contributing to the polarized 
structure functions, 

x 



g{{x) - lfi{x) 



|[ff(s) -■*?(*)] • (28) 



The distributions found in flql , where the nucleon unpolarized structure functions F^x), F£{x) 
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xq(x), and xFz(x) were also described, give the following results 



Au(x) 
Ad(x) 



2.66x 



2.66x 



-0.203 



-0.203 



(1 



(1 



\2.34 



\2.34 



exp 



exp 



1.00 



i -l 



0.235 

x + 0.068 
0.235 



+ 1 



+ 1 



0.123 



exp 



exp 



n -r 



0.235 

x - 0.200 
0.235 



+ 1 



+ 1 



(29) 

in good agreement with the Eq. (28) and with Fit B, as shown in Fig. 5. [] 

We conclude that the SLAC data might be fairly described in terms of Au(x) and Ad(x), with 
a smaller first moment for Au(x) than expected. We cannot, however, conclude that the SLAC 
data give evidence against the Bjorken sum rule. In fact, with r\ u and rjd given in Eqs. (U) and 
(|l7|), it is possible to find a good fit to the distributions by leaving rjo free with a value 1/2 of the 
rhs of Eq. ©. 

Concerning the very important issue of the validity of the Bjorken sum rule, it is worth stressing 
that the SMC data on gf(x) imply a value for I n more negative than the SLAC data with polarized 
He 3 , expecially from the contributions of g™(x) in the small x region not yet measured by SLAC. 
More precise measurements in this region would be very crucial to clarify this point: indeed, two new 
experiments at SLAC, E154 and E155, will be running this year and next year with an electron beam 
energy of 50 GeV/c, and we hope their results will help solving this important problem. Precise 
measurements in the small x region are in fact very relevant also since the negative contribution 
from the gluons advocated to account for the defect in the Ellis and Jaffe sum rule for the proton 
is expected to dominate the small x region. 



1 The values for b u in Eqs. (|l|), (|§) and (||) are smaller than the exponents fixed in Ref. Jl(J by the positivity 
requirement that Au(x) should be smaller than parameterization assumed there for u(x). 

We think there is no problem of positivity since Au(x) as obtained by Eqs. (|^) and ( |23| ) is in good agreement, as 
shown by Fig. 5, with Au(x) = v) (x) — u^(x) of Ref. jl6[ shown in Eq. (p9[), with u(x) — v) (x) + u^x) in good 
agreement with the experimental unpolarized distributions. 
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Figure Captions 



Figure 1. The predictions about xgf(x) corresponding to Fit A (continuous line), Fit B 
(dashed line) and to the values of the parameters reported in Eq. ( |2^ ) (dotted line) are compared 
with SLAC data. 



Figure 2. The predictions about xg™(x) corresponding to Fit A (continuous line), Fit B 
(dashed line) and to the values of the parameters reported in Eq. (|2^ ) (dotted line) are compared 
with SLAC data. 



Figure 3. The predictions about xgf(x) corresponding to Fit A (continuous line), Fit B 
(dashed line) and to the values of the parameters reported in Eq. p^ ) (dotted line) are compared 
with SLAC data. 



Figure 4. The quantities x Au(x) (continuous line), x Ad(x) (dashed line) and x AG(x) (dotted 
line) corresponding to the Fit A are shown. 

Figure 5. The predictions of x Au(x) (continuous line) and x Ad(x) (dashed line) obtained from 
Ref. [16] are compared with x Au(x) (dotted line) and x Ad(x) (dashed-dotted line) corresponding 
to Fit B. 
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TABLE I 





Fit A 


Fit B 


Fit C 


X 2 /Ndf 


1.8 


1.0 


0.63 


a u = a-d 


0.57 ±0.03 


0.98 i 0.09 


0.46 i 0.15 


K 


2.0 ±0.3 


1.8 ±g:i 


3.64 (fixed) 


b d 


3.0 ±0.1 


4.9 i 0.8 


4.64 (fixed) 


7« = Id 


1.0 ±0.8 


0i3 


18.36 i 14.49 




1 (fixed) 




1 (fixed) 


b G 


20 ±1 




7.44 i 3.52 


1G 


Oil 




(fixed) 


Vu/fju 


1 (fixed) 


0.76 i 0.03 


1 (fixed) 


Vd/Vd 


1 (fixed) 


0.93 i 0.08 


1 (fixed) 


vg/vg 


1 (fixed) 


(fixed) 


1 (fixed) 




0.16 


0.17 


0.15 




-0.039 


-0.037 


- 0.042 


AG( 2 ) 


0.10 





0.24 



Table 1: We report, for a comparison, the fitted values for the parameters of the polarized parton 
distributions and their first and second momenta for the cases discussed in the paper. Fit C is the 
one of Ref. [10] based mainly on SMC data. Fit A concerns SLAC data with mainly the same 
ansatz of Ref. [10], while Fit B concerns the same data assuming no gluon anomaly contribution. 
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Fig. 1 




Fig. 2 




Fig. 3 




Fig. 4 




Fig. 5 




